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Abstract 

Allan  Deviation  computations  of  linear  frequency  synthesizer  systems  have  been  reported  pre¬ 
viously  using  real-time  simulations.  Even  though  it  takes  less  time  compared  with  the  actual 
measurement,  it  is  still  very  time  consuming  to  compute  the  Allan  Deviation  for  long  sample  times 
with  the  desired  confidence  level.  Also  noises,  such  as  flicker  phase  noise  and  flicker  frequency 
noise,  can  not  be  simulated  precisely.  The  use  of  frequency  domain  techniques  can  overcome  these 
drawbacks.  In  this  paper  the  system  error  model  of  a  fictitious  linear  frequency  synthesizer  is 
developed  and  its  performance  using  a  Cesium  (Cs)  atomic  frequency  standard  (AFS)  as  a  reference 
is  evaluated  using  frequency  domain  techniques.  For  a  linear  timing  system,  the  power  spectral 
density  at  the  system  output  can  be  computed  with  known  system  transfer  functions  and  known 
power  spectral  densities  from  the  input  noise  sources.  The  resulting  power  spectral  density  can  then 
be  used  to  compute  the  Allan  Variance  at  the  system  output .  Sensitivities  of  the  Allan  Variance  at 
the  system  output  to  each  of  its  independent  input  noises  are  obtained,  and  they  are  valuable  for 
design  trade-off  and  trouble-shooting. 


1.  Introduction 

A  fictitious  linear  frequency  synthesizer  is  used  in  this  analysis.  The  system  generates  a  system 
output  frequency  from  a  reference  Cs  Atomic  Frequency  Standard  (AFS)  operating  at  a  different 
frequency.  The  system  block  diagram  is  shown  in  Figure  1.  A  reference  epoch  is  generated 
every  1.5  s  based  on  the  AFS  frequency  and  another  1.5  s  interval  system  epoch  is  generated 
by  the  system  clock,  a  voltage  control  crystal  oscillator  (VCXO),  Both  epochs  are  input  to  the 
Phase  Meter  (PM),  and  the  PM  computes  the  timing  error  between  them.  Based  on  the  timing 
error  value,  the  loop  adjusts  the  phase  of  the  VCXO-generated  epoch,  so  that  the  VCXO  is 
phase-locked  to  the  reference  AFS. 

In  this  paper  the  system  error  model  of  the  linear  frequency  synthesizer  is  developed  and  the 
performance  at  the  system  output  is  evaluated  using  frequency  domain  techniques.  For  a  linear 
system,  the  power  spectral  density  at  the  system  output  can  be  computed  using  known  system 
transfer  functions  and  known  power  spectral  densities  from  the  three  independent  noise  sources: 
AFS,  VCXO  and  PM.  The  resulting  power  spectral  density  can  then  be  used  to  compute  the 
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Allan  Variance  at  the  system  output.  Compared  with  time  domain  techniques,  the  use  of 
frequency  domain  technique  offers  several  benefits  such  as:  (1)  it  provides  another  independent 
evaluation  of  the  system  performance,  (2)  flicker  noise  can  be  implemented  precisely,  (3) 
computation  time  for  the  Allan  Variance  of  a  given  sample  time  (r)  is  very  short  and  is  roughly 
the  same  for  either  short  or  long  sample  time  (r).  The  last  two  items  are  major  drawbacks  for 
the  real-time/Monte  Carlo  simulations.  Sensitivities  of  the  output  Allan  Deviation  to  each  of 
its  input  noises  are  provided,  and  they  are  valuable  for  design  trade-off  and  trouble-shooting. 

2.  Noise  Models 

2,1  Phase  Meter  Noise 

The  PM  noise  with  resolution  of  100  ps  is  modeled  as  a  white  process  with  a  constant  power 
spectral  density  for  all  frequencies 


Vp(f)  = 


(1  x  10~10)2 
6 


And  the  Allan  Variance  of  the  PM  noise  is^l: 


tp(t)  = 


O'!  2 


(1  x  10]U) 
6t2 


(1) 


(2) 


2.2  AFS  Noise 

The  Allan  Variances  for  the  AFS  is  specified  as: 


<%(*) 


9.0  x  10~22 
r 


+  1.0  x  102° 


(3) 


This  is  consistent  with  that  of  the  commercially  available  HP  5071  Cs  AFS. 

We  assume  that  the  two  terms  are  independent  of  each  other.  Thus  we  can  compute  the 
associated  noise  power  spectral  densities  using  well  known  techniquestO: 


SA(f)  =  1.8  x  10-21 


+ 


7.2134  x  10-27 

7 


(4) 


2.3  VCXO 

The  VCXO  frequency  noise  is  assumed  to  have  the  following  Allan  Variance: 

ct2(t)  =  10- 24  +  10~27t  .  (5) 
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Similarly,  we  assume  that  the  two  terms  of  equation  (5)  are  independent  of  each  other.  As 
above,  the  power  spectral  density  of  the  VCXO  frequency  noise  can  expressed  asPl; 


SJf)  = 


7.2134  x  lO-25  1.519  x  10-28 

1 _ 


(6) 


3.  System  Transfer  Functions 

The  system  model  as  shown  in  Figure  1  results  in  the  equivalent  system  error  model  indicated 
in  Fig.  2  obtained  using  the  Z-transform  formalism.  Ts  is  the  sample  period  of  the  system, 
and  is  1.5  s.  Also  two  delays  of  one  epoch  each  are  introduced  in  Figure  2  to  account  for  the 
fact  that  the  effects  of  the  computed  VCXO  frequency  modification  in  the  current  epoch  will 
not  show  on  the  Phase  Meter  until  two  epochs  later.  The  transfer  functions  relating  the  system 
output  frequency  (T)  to  the  input  noises  of  the  AFS(A),  VCXO(V)  and  Phase  Meter  (P)  are: 


SYSTEM  OUTPUT 

PM 

11 

3 

ii 

T{Z) 

P(Z) 

U) 

SYSTEM  OUTPUT 

AFS 

=  Hta(Z)  = 

T(Z) 

A(Z ) 

U) 

SYSTEMOUTPUT 

=  Htv(Z)  = 

T{Z) 

(0) 

VCXO 

V(Z) 

The  Bode  Plots  (frequency  responses)  for  these  second  order  transfer  functions  with  a  time 
constant  of  50  s  are  shown  in  Figure  3.  It  is  seen  that  Wta(Z)  is  a  lowpass  filter,  H tv(Z)  is  a 
highpass  filter,  and  Hrp(Z)  is  a  low  gain  highpass  filter  to  reduce  the  PM  quantization  noise. 

4.  Power  Spectral  Density  of  the  System  Output  Frequency 

The  power  spectral  density  of  the  system  output  (ST(f))  can  be  computed  asl2l: 


St(I)  —  StaU )  +  Srv(f)  +  (10) 

where: 

StaU)  =  \HTA(U2*fTil)\2  *  SA(f),  due  to  the  AFS  noise  (11) 

StvU)  ~  \HTv{U2lTfTs)\2  x  Sv(f ),  due  to  the  VCXO  noise  (12) 


and 
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Stf(I)  =  x  5p(/),  due  to  the  PM  noise. 


(13) 


The  power  spectral  densities  for  the  noises  (SaU),  Sv(f)  and  Sp(f))  and  their  contributions 
to  the  system  frequency  output  (Sta(I),  Stv(/)  and  Spp(f))  were  calculated  using  equations 
(1),  (4),  (6),  (11),  (12)  and  (13),  and  are  shown  in  Figure  4.  As  can  be  seen  the  system  output 
power  spectral  densities  are  shaped  by  their  corresponding  transfer  functions.  Even  though  the 
power  spectral  densities  of  the  PM  ( Sp(f ))  and  the  VCXO  (Sv(/))  are  large  when  compared 
with  that  of  the  AFS,  their  contributions  to  the  system  output  ( SppiJ )  and  STv(f))  have  been 
greatly  reduced  by  their  corresponding  transfer  functions  (Hrp{j2nfTs)  and  Hpv 0"2tt /Ts)) 
especially  at  low  frequency.  Similarly,  contributions  from  the  AFS  (Sj^)  is  also  suppressed 
considerably  for  high  frequencies. 


5.  Allan  Deviation  of  the  System  Output  Frequency 

The  temporal  behavior  of  the  timing  or  frequency  system  is  normally  characterized  by  the  Allan 
Deviation  or  square  root  of  the  Allan  Variance,  <t2(t),  where  r  is  the  frequency  sample  time. 
The  Allan  Variance  is  related  to  the  power  spectral  density  by  [3] 


a2y(r)  =  2j^Sy(f)AD(fT)df 

where: 


(14) 


Ad(Ft) 


sin4(7r  fr) 

(*  M2 


(15) 


and  fn  is  the  Nyquist  frequency  of  the  system  and  is  equal  to  1/3  Hz  for  a  sampling  period 
(T,)  of  1.5  s. 


5.1  Computation  Consideration 

In  many  cases  integration  of  the  equation  (14)  can  not  be  carried  out  analytically,  so  it  must  be 
done  by  computer  using  numerical  integration.  It  is  important  to  choose  a  proper  integration 
step  size  to  achieve  the  desired  accuracy  in  the  computation. 

The  magnitude  of  the  oscillatory  window  function  (AD(/))  as  given  in  equation  (15)  is  inversely 
proportional  to  the  square  of  /  for  a  given  r,  and  is  plotted  in  Figure  5  for  r  =  100  s  and 
r=1000  s.  As  can  be  seen,  its  magnitude  decreases  rapidly  after  a  few  periods  and  its  bandwidth 
decreases  as  t  increases.  Figure  5  shows  that  the  bandwidth  is  roughly  10~s  Hz  for  r=1000  s. 
and  10-2  Hz  for  r- 100  s.  The  portion  of  the  power  spectral  density  of  the  system  noise  outside 
the  bandwidth  of  the  window  function  has  negligible  effect  of  its  Allan  deviation  computed  by 
using  equation  (14). 
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Given  the  limitations  in  available  computer  memory  it  was  found  by  trial  and  error  that  to 
provide  adequate  numerical  accuracy,  integration  of  the  equation  (14)  can  be  carried  out  by 
using  20  integration  steps,  either  for  the  first  50  periods  of  the  window  function  or  up  to  the 
Nyquist  frequency,  whichever  is  smaller.  This  technique  is  valid  for  sample  times  of  up  to 
100,000  s  provided  that  the  power  spectral  density  function,  Sy(f),  does  not  contain  any  white 
phase  modulation  noise  component,  whose  frequency  power  spectral  density  is  proportional  to 
f2.  Fortunately  this  condition  is  met  by  the  frequency  power  spectral  density  at  the  system 
output.  For  white  phase  modulation  noise  such  as  the  PM  noise,  the  upper  limit  of  the 
numerical  integration  has  to  be  set  at  the  Nyquist  frequency.  For  larger  sample  times,  r  > 
100,000  s,  a  smaller  integration  step  size  is  needed.  The  average  time  to  compute  the  Allan 
Deviation  for  a  given  jtau  is  less  than  15  s  using  a  PC  with  Intel  486  DX2/50  CPU. 

To  show  that  this  numerical  integration  technique  is  accurate  enough  the  Allan  Deviations  of 
the  two  noise  sources  (AFS  and  VCXO)  are  computed  from  their  corresponding  power  spectral 
densities  using  equation  (14)  and  are  shown  in  Figure  6.  Figure  6  also  shows  the  corresponding 
specified  Allan  Deviations,  from  which  the  power  spectral  densities  were  derived,  as  discussed 
in  section  2.  As  can  be  expected,  the  computations  are  in  very  good  agreement  with  the 
respective  specifications.  The  Allan  Deviation  of  the  PM  noise  is  also  depicted.  Figure  6  shows 
that  the  PM  Allan  Deviation  is  predominant  for  sample  times  up  to  30  s  and  that  of  the  VCXO 
noise  predominates  for  longer  time. 

5.2  System  Allan  Deviation 

The  contribution  to  the  Allan  Deviations  at  the  system  output  (crjviCr),  <jtp(t)  and  ctTv(t)) 
due  to  each  of  the  independent  input  noises  are  computed  using  equation  (14)  with  their  power 
spectral  densities  at  the  system  output  Stp(I),  and  Stv{I))-  The  results  are  shown 

in  Figure  7  and  they  can  be  considered  as  the  sensitivities  of  the  system  output  for  each  of  the 
input  noises.  This  technique  can  be  used  very  effectively  during  the  design,  development  and 
testing  phase  of  the  system  to  determine  the  loop  time-constant,  to  define  noise  specifications 
and  to  provide  data  for  trouble-shooting,  ft  will  be  used  below  to  identify  the  causes  of  some 
exceeding  reference  AFS  specification  conditions.  It  is  apparent  from  Figure  8  that  for  short 
sample  time  the  system  performance  is  dominated  by  the  PM,  while  for  long  sample  time  the 
system  performance  is  governed  by  the  AFS.  The  crossover  sample  time  is  around  20  s.  The 
Allan  Variance  of  the  resulting  system  output  can  be  obtained  as: 


ar(T)  =  <7ta(t)  +  <4v(r)  +  v'tp(t)-  (16) 

The  Allan  Deviations  of  the  system  output  (<jt(t))  and  the  reference  AFS  specification  are 
plotted  in  Figures  8.  Figure  8  shows  that  the  Allan  Deviation  of  the  system  output  barely 
exceeds  the  reference  AFS  specification  for  sample  time  from  70  s  to  1000  s.  By  examining 
the  system  sensitivities  as  provided  in  Figure  7,  it  is  found  that  this  condition  is  caused  by  the 
AFS  noise. 
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6.  System  Performance  Using  Other  System  Configurations 

In  the  previous  sections  we  only  use  Tc  ~  50  s,  Allan  Deviation  results  using  other  time 
constants  can  be  obtained  easily.  The  result  for  a  time  constants  of  15  s  is  shown  in  Figures  9. 
Figure  9  indicates  that  the  system  exceeds  the  reference  AFS  Allan  Deviation  specification  for 
r  1000  s.  The  cause  can  be  determined  by  examining  the  sensitivities  as  shown  in  Figure  10, 
and  is  identified  to  be  the  PM  noise,  which  is  the  dominant  noise  contribution  for  short  x;  its 
effect  at  the  system  output  is  not  suppressed  enough  for  the  system  with  a  short  time  constant. 
To  reduce  the  effect  due  to  the  PM  for  short  sample  time,  a  better  PM  with  resolution  of  20 
ps  was  used  and  the  result  is  shown  in  Figure  11.  As  can  be  seen  the  Allan  Deviation  at  the 
system  output  is  less  than  that  of  the  reference  AFS  for  short  sample  time. 

7.  Conclusion 

An  efficient  method  of  computing  the  Allan  Deviation  at  the  output  of  a  linear  system  with 
known  input  power  spectral  densities  is  presented.  Since  the  computation  time  is  not  a  function 
of  the  sample  time  (fau),  this  technique  is  very  attractive  to  compute  Allan  Deviations  for 
long  sample  times.  Sensitivities  of  the  Allan  Deviation  at  the  system  output  for  each  of  its 
independent  input  noises  are  also  provided  and  they  are  valuable  for  design  trade-off  and 
trouble-shooting.  Potential  situations  in  which  the  system  could  exceed  the  reference  AFS 
specification  are  pointed  out,  and  causes  are  identified. 
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Figure  1.  Linear  Frequency  Synthesizer 


Figure  2.  System  Error  Model 
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Figure  5.  Plot  of  Window  Function  as  functions  of  f  and  tau 


Figure  6.  Allan  Deviations  of  PM,  VCXO  and  AFS 
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